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PLAN DU COURS 


OBJECTIF: Conception des algorithmes corrects et efficaces 

PLAN 

□ RAPPELS : NOTATIONS ALGORITHMIQUES 

□ COMPLEXITE 

□ ALGORITHMES ITERATIFS DE TRIS 

□ RECURSIVITE 

□ DIVISER POUR RESOUDRE 

□ PREUVE D’ALGORITHMES 
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Notations algorithmiques 


□ Un type est un ensemble de valeurs sur lesquelles on definit 
des operations. 

• Types de base : 

y Entier : Operateurs arithmetiques +, *, div, mod 

y Reel : Operateurs arithmetiques +, *, / 

y Booleen : Operateurs logiques et, ou, non 

y Caractere : constante (lettre imprimable) entre 
apostrophe. 

_ Les operateurs relationnels permettant de faire des 
comparaisons: < f < f =, >, >, ^ 

Le Resultat de la comparaison est une valeur booleenne. 
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Notations algorithmiques 


□ Une variable possede : 

- un nom 

- une valeur 

- un type 

(la valeur d’une variable peut changer au cours de 
I’execution) 

Declaration : <variable> : <type> 


□ Une expression, pour un type, est soit une constante, 
soit une variable, soit constitute a I’aide de constantes, 
de variables, de parentheses et des operateurs 
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Notations algorithmiques 
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□ INSTRUCTIONS 

■ Affectation : <variable> := <expression> 


■ Condition : 


si <condition> alors 
action 


fsi 


OU : 


si <condition> alors 

action 1 


sinon 


action 2 


fsi 

( condition est une expression a valeur booleenne; 
action est une instruction ou un bloc destructions separees par ;) 
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Notations algorithmiques 


- Iteration 

• boucle Pour 


pour <variable> := <initiale> a <final> faire 
action 

fpour 


(Ou initiale et finale sont des expressions de meme type que celui de 
la variable controlant la boucle, le type peut etre entier, caractere 
ou enumere) 


Remarque: la boucle pour affecte la valeur de initiale a variable et 
compare cette valeur a celle de finale avant d’executer action. 
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□ Exemple: calcul de 1 +2 + ...+n (n entier> 1 
fixe) 

Programme somme_des_n_premier$Termes 
// partie declaration 

n : entier; s : entier; i : entier; (ou n, i, s : entier) 
Debut 

/ / Lecture des donnees 

Ecrire (" n = ? " ); lire(n); 

/ / calcul de la somme 
s := 0; 

pour i := 1 an faire 
s := s + i; 

fpour; 

/ / affichage du resultat 

ecriref'l + 2 + ... + n = ", s); 


fin 


Notations algorithmiques 


en C 


main ( ) { 

// declaration des variables 
int i, n, s; 

/ / lecture des donnees 

prinf(" n = ? "); scanf(" %d" , &n); 
// calcul de la somme 
s = 0; 

for(i — 1 ; i <— n; i++) s = s + i; 

/ / affichage 

printf(" 1 +2 + ...+ %d = %d \n" , n, s) 

} 
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Notations algorithmiques 


- Iteration 

• boucle tantque • boucle repeter 


tantque < condition> faire 
action 

ftantque 


repeter 

action 

jusque < condition> 


Remarques: 

Condition est une expression a valeur booleenne, cette expression 
doit etre modifiee dans le corps de la boucle (dans action). 

La boucle pour peut etre traduite en boucle tantque (ou en repeter) 
mais I’inverse n’est pas toujours vrai. 
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Programme somme_des_n_premiersTermes 
// partie declaration 

n : entier; s : entier; i : entier; (ou n, i, s : entier) 
Debut 

/ / Lecture des donnees 

Ecrire (" n = ? " ); lire(n); 

// calcul de la somme 
s := 0; 

I:= 1; 

tantque i <= n faire 
s := s + i; 
i := i + 1; 
ftantque; 

// affichage du resultat 

ecriref'l + 2 + ... + n = ", s); 
fin 


Notations algorithmiques 


□ en C 
main ( ) { 

// declaration des variables 
int i, n, s; 

/ / lecture des donnees 

prinf(" n = ? "); scanf(" %d" , &n); 
// calcul de la somme 
s = 0; 

i = 1; 

while (i <= n) { 

s = s + i; 
i = i + 1; 

} 

/ / affichage 

printff' 1 +2 + ...+ %d = %d \n" , n, s) 

} 


FSR Dept.lnfo. 1 5 



Notations algorithmiques 


□ Algorithme 

■ Description formelle d’un procede de calcul qui permet, 
a partir d’un ensemble de donnees, d’obtenir des 
resultats. 

• Succession finie d’operations 

DONNEES algorithme, RESULTATS 


Cheminement d suivre : 


Probleme 


analyse 


codification 

_ Alaorithme _ programme 
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Notations algorithmiques 


D Conception structuree des algorithmes 

Decouper I’algorithme en sous algorithmes plus 
simples. Chaque sous algorithme est un algorithme. 

> Fonctions 

Une fonction est un sous algorithme (sous programme) qui, a partir 
de donnees, retourne un seul type de resultat . 

Une fonction 

• possede un nom 

• communique avec I’exterieur par le biais des parametres 

• retourne un resultat par I’instruction retourner(expression) 
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Notations algorithmiques 


□ Schema d’une fonction : 

fonction <nom de la fonction>(Liste des parametres) : <type du resultat> 
// declaration des variables locales 

debut 

/ / corps de la fonction qui contient I’instruction retourner 

fin 


Les parametres de la definition d’une fonction (appeles formels) sont : 

• types 

• separes par 7 s’il y en a plusieurs 

Les parametres formels sont utilises pour ne pas lire les donnees dans une 
fonction. 

Les variables declarees dans une fonctions (y compris les parametres formels) 
sont appelees variables locales. 
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□ Exemple: factoriel n en 
pseudo code: 

Fonction fact(n : entier) : entier 
m, i : entier; 
debut 

m := 1 ; 

pour i :—2 a n faire 
m := m * i; 

fpour 

retourner(m); 


fin 


Notations algorithmiques 


□ n! en C : 

unsigned int fact (unsigned int n) { 
unsigned int i, m; 
m = 1 ; 

for (i = 2; i<= n; i++) 
m = m * i; 
return m; 

} 
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Notations algorithmiques 
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o L’appel d’une fonction est utilise dans une instruction sous la forme : 

<nom de la fonction> (liste des parametres effectifs) 

> Les parametres formels et effectifs doivent correspondre en nombre 
et en type 

> Lors d’un appel : 

• les parametres formels regoivent les valeurs des parametres 
effectifs correspondant 

• le corps de la fonction est execute jusqu’au premier retourner 
rencontre 

• I’execution se poursuit (a I’dresse de retour) dans la fonction 
appelante. 

Remarque: le type de retour d’une fonction peut etre vide, et dans 
ce cas on ecrit retournerQ ou pas d’instruction retourner; 
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□ Exemple: calcul de 1 +2 + ...+n (n 
entier> 1 fixe) 

Programme somme_des_n_premiersTermes 
// partie declaration 
n : entier; r : entier; (ou n, r : entier) 
Debut 

/ / Lecture des donnees 

Ecrire (" n = ? " ); lire(n); 

/ / appel et utilisation de la fonction 
sommeArith 

r := sommeArith(n); 

/ / affichage du resultat 

ecriref'l + 2 + ... + n = ", r); 


fin 


Notations algorithmiques 


□ Fonction pour calculer la somme 

1 +2 + ...+m 



Fonction sommeArith ( m : entier) : entier 
i, s : entier; 
debut 
s := 0; 

pour i := 1 dm faire 

S := S + i; 

fpour; 

retourner(s); 

fin 
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Notations algorithmiques 


o Deux types de passages des parametres : 

• Passage par valeur: la fonction travaille sur une copie du 
parametre effectif transmis; i.e. la valeur du parametre effectif 
n’est pas modifiee apres I’appel de la fonction. 


• Passage par adresse (ou par reference): 

- I’identificateur du parametre formel est precede par le mot ref. 

- la fonction travaille directement sur I’identificateur du parametre 
effectif; i.e. toute modification sur le formel entraine la meme 
modification sur le parametre effectif correspondant. 
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Passage par valeur 

Fonction echnger(x: reel, y : reel) : vide 

z : reel; 

debut 

Z := x; 

X := y; 

y == z; 

fin 

Fonction appelanteQ 

a, b : reel; 

debut 

a := 2; b:= 7; 
echanger(a,b); 
ecrire( a = , a); 
ecrire( b = , b); 
fin 


Les resultas affiches par la fonction appelante : 

a = 2 b = 7 


Notations algorithmiques 


Passage par reference 

Fonction echnger(ref x: reel, ref y : reel) : vide 

z : reel; 

debut 

Z := x; 

X := y; 

y == z; 

fin 

Fonction appelanteQ 

a, b : reel; 

debut 

a := 2; b:= 7; 
echanger(a,b); 
ecrire( a = , a); 
ecrire( b = , b); 
fin 

Les resultas affiches par la fonction appelante : 

a = 7 b = 2 
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Notations algorithmiques 



□ Un alqorithme se comporte comme une fonction sauf 
que Ton ne s’occupe pas des declarations des variables 
ni de leurs types. 


> Schema d’un alqorithme : 

Algorithme <nom de l’algorithme> 

Donnees : // 1 es variables qui sont des donnees de I’algo. 

<nom de l’algorithme>(Liste des parametres) Resultat(s): // variable(s) contenant le(s) resultat(s) 


debut 

/ /bloc d’instructions Ou 

fin 


debut 

/ /bloc d’instructions 

//ne contenant pas retourner 

fin 


On omet la partie declaration des variables locales en adoptant la 
regie: les variables simples sont en minuscule et les tableaux en 
majuscule. 
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Notations algorithmiques 


□ Algorithme pour calculer n! 

factoriel(n) 

// n > 0 

debut 

m := 1 ; i := 1 ; 
tantque i < n faire 
i := i + 1; 
m := m * i; 
ftantque 
retourner(m); 
fin 
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TABLEAUX STATIQUES 


□ TABLEAUX 

■ Un tableau est utilise pour representer une suite d’elements de meme type 
( T — (t] /^ 2 '* * v^ n ))* 

> Declaration d’un tableau (a un dimension): 

<nom du tableau> : tableau [1 .. max] de <type des elements> 

(exemple T : tableau [1 ..20] de reel) 

ou : - max est une constante entiere (positive) 

- le type des elements est quelconque (de base ou declare). 

• La taille (ou longueur) d’un tableau est le nombre d’elements du tableau 
elle est toujours inferieure ou egale a max. 

• les tableaux, en algorithmique, commencent a I’indice 1 (en C et en java, 
ils commencent a I’indice 0) 
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TABLEAUX STATIQUES 


> Operations (de base) sur les tableaux: 

❖ Acces a un element du tableau : 

<nom du tableau> [<indice>] 

(indice est une expression de type entier) 

❖ Parcours : (on utilise un indice et une boucle pour ou tant que) 
(exemple : 

pour i= Ian f aire //n est le nombre d’elements du tableau T 

traiter( T[i] ) // traiter() est une fonction ou algorithme a definir 

fpour; 

❖ Recherche d’un element dans un tableau 

❖ Insertion d’un element 

❖ Suppression d’un element 
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TABLEAUX STATIQUES 
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Exemple : recherche de la position d’un element dans un tableau de reels. 


fonction localiser(T: tableau[l ..max]de reel, n : entier, val :reel) : entier 
i : entier; 
trouve : booleen; 
debut 

i := 1 ; trouve := faux; 
tantque (i < n) et (non trouve) faire 
si (T[i] = val) alors 
trouve := vrai; 
sinon i := i + 1 ; 
fsi; 

ftantque; 

si (trouve) alors retourner(i) 
sinon retourner(O); 
fsi 
fin 
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Algo.ll 


□ insertion d’element x dans un tableau Tan elements a la position p. 

Fonction insererfT : tableau[l ..max] de reel, ref n : entier, x : reel, p : entier) : vide 
//I <p < n 
i : entier; 
debut 
i := n; 

tantque i > p faire 

T[i+1]:= T[i]; 

i := i — 1 ; 
ftantque 

T[p] := x; 

n := n + 1 ; 
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Algo.ll 


Suppression d’une valeur d’un tableau T, qui se trouve a la position 
P- 

Fonction supprimer(T : tableau[l .. max] de reel, ref n : entier, p : entier) : vide 
i : entier; 
debut 
i := p; 

tantque i < n faire 
T[i] := T[i+ 1 ]; 

i := i + 1 ; 
ftantque; 
n := n - 1 ; 
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Algo.ll 


o Remarques : 

^ Le nom du tableau, dans une liste de parametres 
formels, est une reference. 

(en C, I’adresse de T[i], notee T + i, est calculee 

adresse(T[i]) = adresse(T[0]) + sizeof(<type des elements>) * i) 

^ L’insertion (resp. suppression) d’un element a un 
indice i, necessite un decalage des (n- i +1 ) 
derniers elements d’une position a droite (resp. a 
gauche) 
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STRUCTURE (ENREGISTREMENT) 

□ Structures (ou enregistrements) 

■ Le type structure est utilise pour representer une suite d’elements pas 
necessairement de meme type, chaque element est appele champs. 

> Declaration d’un type structure: 

<nom de type structure> = structure 

<variable_champsl > : <type_champsl >; 

<variable_champs2> : <type_champs2>; 

fstructure 

> Declaration d’une variable de type structure: 

<variable_structure> : <nom de type structure> 
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STRUCTURE (ENREGISTREMENT) 


Par analogie au type struct de C: 

Declaration : 

struct <nom de la structure> { 

<type_champsl > <variable_champsl > 
<type_champs2> <variable_champs2> 


} 

Declaration d’une variable de type structure: 
struct <nom de la structure> <variable> 

Ou avec la definition de type : typedef 

typedef struct <nom de la structure> <type_structure>; 

Declaration d’une variable : 

<type_structure> <variable> 
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STRUCTURE (ENREGISTREMENT) 



□ Exemples : 

Adresse = structure 

numero_rue : entier; 
nom_rue : tableau[l ..20] de caractere; 
code_postal : entier; 
fstructure 

Point = structure 
abscisse : reel; 
ordonnee : reel; 
fstructure 


□ En C : 

struct adresse { 

int numero_rue; 
char[20] nom_rue; 
int code_postal; 

} 

typedef struct adresse Adresse; 

typedef struct point { 

float abscisse, ordonnee; 

}Point; 


Declaration des variables: 
adr : Adresse; 
p : Point 


/ /Declaration des variables: 
Adresse adr; 

Point p; 
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STRUCTURE (ENREGISTREMENT) 

> Operation sur les structures : 

❖ Acces a un champs : 

<variable de type structured <variable_champs> 
(ex: p.abcisse designe I’abscisse du poin p) 

❖ Affectation : 

<variable_type_structure> := <variable_type_structure> 

(L’affectation se fait champs par champs) 

Exemple : manipulation des complexes 
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STRUCTURE (ENREGISTREMENT) 


■ Declaration de type complexe: 

Complexe = structure 

p_reel / p_imag : reel; 

fstructure 

fonction plus(z1 : Complexe, z2 : Complexe) : Complexe 

z : Complexe; 
debut 

z.p_reel := zl.p_reel + z2.p_reel; 
z.p_imag := zl.p_imag + z2.p_imag; 
retourner(z); 
fin 


□ Declaration de type complexe en C: 
typedef struct complexe { 
double p_reel, p_imag ; 

} Complexe; 

Complexe plus(Complexe zl, Complexe z2){ 
Complexe z; 

z.p_reel = zl.p_reel + z2.p_reel; 
z.p_imag = zl.p_imag + z2.p_imag; 
retun z; 

} 
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ALGORITHMIQUE II 


NOTION DE COMPLEXITE 



SMI AlgoII 



Complexity 


□ Comment choisir entre differents algorithmes pour 
resoudre un meme probleme? 

□ Plusieurs criteres de choix : 

> Exactitude 

> Simplicity 

> Efficacite (but de ce chapitre) 
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Complexity 


□ [.’evaluation de la complexity d’un algorithme se fait 
par I’analyse relative a deux ressources de I’ordinateur: 

o Le temps de calcul 

o L’espace memoire, utilise par un programme, pour 
transformer les donnees du probleme en un ensemble 
de resultats. 

L’analyse de la complexity consiste a mesurer ces deux 
grandeurs pour choisir I’algorithme le mieux adapte 
pour resoudre un probleme.(le plus rapide, le moins 
gourment en place memoire) 


On ne s’interesse, ici, qu’a la complexity temporelle c.a d. qu’au temps de calcul 
(par opposition a la complexity spatiale) SMI A| g o11 




Complexity 


□ Le temps d’execution est difficile a prevoir, il peut 
etre affecte par plusieurs facteurs: 

> la machine 

> la traduction (interpretation, compilation) 

> I’environnement (portage ou non) 

> L’habilete du programmeur 

> Structures de donnees utilisees 
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Complexity 
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□ Pour pallier a ces problemes, une notion de 
complexity plus simple, mais efficace, a ete definie 
pour un modele de machine . Elle consiste a 
compter les instructions de base executees par 
I’algorithme. Elle est exprimee en fonction de la 
taille du probleme a resoudre. 

> Une instruction de base (ou elementaire) est soit 
une affectation, un test, une addition, une 
multiplication, modulo, ou partie entiere. 
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Complexity 


□ La complexity depend de la taille des donnees de I’algorithme. 

■ Exemples : 

> Recherche d’une valeur dans un tableau 

— > taille (= nombre d’elements) du tableau) 

> Produit de deux matrices 
— > dimension des matrices 

> Recherche d’un mot dans un texte 

— > longueur du mot et celle du texte 
On note generalement: 

n la taille de donnees, T(n) le temps (ou le cout) de I’algorithme. 
T(n) est une fonction de IN dans IR + 
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Complexity 


□ Dans certains cas, la complexity ne depend pas 
seulement de la taille de la donnee du probleme 
mais aussi de la donnee elle-meme. 

■^Toutes les donnees de meme taille ne generent 
pas necessairement le meme temps d'execution. 

— > (Ex. la recherche d’une valeur dans un tableau 
depend de la position de cette valeur dans le 
tableau) 
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Complexite 


□ Une donnee particuliere d’un algorithme est appelee instance 
du probleme. 

□ On distingue trois mesures de complexite: 

1 . Complexite dans le meilleur cas 

T Min (n) — min {T (d) ; d une donnee de taille n } 

2. Complexite dans le pire cas 

T Mx (n) — max {T(cf) ; d une donnee de taille n } 

3. d ans la cas moyen 

T MOY (n) “ 2] pfcf).T(cf) 

d de taille n 

p(d) : probability d'avoir la donnee d 
T MIN (") ^ T MOY ( n > ^ T MAx( n > 
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□ Exemple. Complexite de la 
recherche d’un element x dans un 
tableau A a n valeurs. 


i := i 

tantque (i <n) et (A[i] ^x) faire 
i := i+i; 

jbmtque 

si i>n alors retoumer(faux ) ; 
sinon retoumer(vrai); 

fa 

On note par: 

a : le cout d’une affectation 
t : cout d’un test 
d : cout d’une addition 


Complexite 


□ Cas le plus favorable , x est le 

premier element du tableau: 

T m» = 1 a + 3t 

□ Pire des cas . x n’est pas dans le 
tableau: 

T max( n ) - (n+1 )a +(2n+2)t +nd 

□ En movenne : 
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Complexite 


■ Complexite en moyenne : 

□ On note par : 

D, (1 <i< n) : ensemble de donnees (de faille n) ou x est 
present a la i eme position 

D n+1 : ensemble de donnees ou x n’est pas present 

□ On suppose que la probability de presence de x dans une 
donnee est q. De plus, dans le cas ou x est present, on 
suppose que sa probability de presence dans I’une des 
positions est de 1 /n 

On a : 

P(Di) = q/n, T(Dj) = i a +(2i+l)t + (i-l)d ; (1 <i< n) 
p( D n+1 ) = 1 - q , T(D n+1 )=T max = (n+1 )a+(2n+2)t+nd 
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Complexity 



p(Dj) = q/n , T(D;) = i a +(2i+l )t + (i-1 )d ; (1 <i< n) 

P( D n+1 ) — 1 - q, T(D n+1 )— T max (n) = (n+1 )a+(2n+2)t+nd 


T MOY (n) “ 2 p(cf).T(cf) 

d de faille n 

p[d) : probability d'avoir la donnee d 

n 

Tmoyin] = y p(Di)T{Di) + p(Dn + i)T(Pn + 1) 

£=1 


T* 


Tmoyin) = — / (ia + (2i + l)t + (i — 1) d) + (l — q)[(?i 4- 1 )a + (2 n + 2 )t + nd] 
nZ_ i 


E — 1 


VI 4 1 71 — 1 

Trnoy'ji) — q [ — - — a 4- (n 4 2)t H — id] 4 (l — q) [(n 4- l) a 4- (2n 4- 2)t + r.d] 


Cas ou q=l , i.e. x est toujours presenttf’moyOrO = l / 2 : [(n + 1 )o + (2 n + 4}t + (n - l)rf] 

"sMjriigo# 2tr+dJ. + V 2 (° + -d) 


On remarque que la complexity de cet algo, est 
de la forme: a n 4 p , a et P son des constates 
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Complexite asymptotique 
Comportement de T(n) 


□ Pour mesurer la complexite d’un algorithme, il ne s’agit pas de faire un decompte exact du 
nombre d’operations 1(11), mais plutot de donner un ordre de grandeur de ce nombre pour n 
assez grand. 

□ Notation de Landau 

❖ "grand O" 

T(n) = 0(f(n)) ssi 

3 c > 0 3 n 0 > 0 Vn>n 0 T(n)< c.f(n) 

❖ "grand omega" 

T(n) = £1 ( f(n )) ssi 

3 c > 0 3 n 0 > 0 Vn>n 0 T(n)> c.f(n) 

❖ "grand theta" 

T(n) = 0(f(n)) ssi 

3c ? >0 3c 2 >03n 0 >OVn>n 0 c n f (n)< T(n)< c 2 f(n) 


Remarque: les constantes c, c u c 2 et n 0 sont independantes den 
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Complexity 



D’une maniere generale, f : IR — > IR 

f(x) = O(g)x)) s’il existe un voisinage V de x 0 et une constante k>0 tels 
que |f(x)|< k|g(x)| , (x e V) 


- Si la fonction g ne s’annule pas, il revient au meme de dire que le 

. /'CO 

rapport 


£f(0 


est borne pour x e V. 
-Exemple: au voisinage de 0, on a: 
x 2 = O(x), ln(l+x) = O(x) 


- Au voisinage de I’infini (comme pour le cas de la complexity), il 
existe a>0 (V = ]a, + [) et k>0 t.q 

|f(x)|< k|g(x)| , Vx > a 


et on dit que f est dominee asymptotiquement par g. 
(Au voisinage de +°° , on a: x = 0(x 2 ), In x = O(x)) 
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Complexity 



□ Remarques 


1. 

2 . 

3 . 

4 . 


f = O(g) <=> 
de I’infini. 


le quotient 


est bornee au voisinage 


II im 

SP — +C-Z 


fix) I 

&o^y\ 


a — * f = Ofej 


si a — 1, on ecrit f ~ g et on a f — 0(g). 

f = 0(g) ne signifie pas que le quotient f(x)/g(x) 
tend vers une limite (1 en particulier). 

Exemple. f(x) = x(2+sin x). On a x < f(x) < 3x 
Vx>0, done f(x) = 0(x). En revanche, le quotient 
f(x)/x ne tend vers aucune limite lorsque x — » +°o 
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Complexity 


□ Abus de notation: 

On a par definition: 

O(f) = {g: IN— HR/3c>0 3 r> 0 > 0 g(n) < c.f(n) , V n > n 0 } 

H(f) = {g: IN — HR / 3 c > 0 3n 0 >0 g (n) > c.f(n) , V n > n 0 } 

0 (f) = O(f) n n(f) 

La difficulty, dans la familiarisation avec ces concepts, provient de la convention 
de notation (de Landau) qui veut que Ton ecrive : 

g = O(f), ou encore g(n) = 0(f(n)) au lieu de g e O(f) 

De maniere analogue, on ecrit O(f) = O(g) lorsque 0(f)cz O(g) 

(il en est de meme pour les notations £1 ou 0) 
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Complexity 


□ Exemple. Soit la fonction T(n) = Vi n 2 + 3n 

> T(n) = £2(n) ( n 0 = 1 , c = Vi) 

> T(n) = 0(n 2 ) (n 0 = 1, q = Vi, c 2 = 4) 

> T(n) = 0(n 3 ) (n 0 = 1, c = 4) 

> T(n) ^ O(n) 

Supposons que T(n) = O(n) 

3 c >0, 3n 0 >0 : V 2 n 2 +3n < cn \/n >n 0 
done c > V 2 n , contradiction, (la constante c ne peut dependre de n) 
Remarques 

Si T(n) est un polynome de degre k alors T(n) = 0(n k ) 

2. 0(n k ) d 0(n k+1 ) V K > 0 (idem pour 0) 

0(f(n)) cz 0(f(n)) pour toute fonction f positive 
0(1 ) utilise pour signifier « en temps constant » 
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Complexity 


□ Remarques pratiques: 

Le cas le plus defavorable est souvent utilise pour 
analyser un algorithme. 

La notation O donne une borne superieure de la 
complexity pour toutes les donnees de meme 
taille(suffisamment grande). Elle est utilisee pour evaluer un 
algorithme dans le cas le plus defavorable. 

T(n) < cf(n) signifie que le nombre d’operations ne peut 
depasser cf(n) iterations, pour n’importe quelle donnee de 
longueur n. 

Pour evaluer la complexity d’un algorithme, on cherche un 
majorant du nombre d’operations les plus dominantes. 

Dans les notations asymptotiques, on ignore les constantes. 

L’alqorithme de recherche dans un tableau a n elements, cite precedemment. est en O(n) 
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Complexity 


□ Ordre de grandeur courant 

0(1) : complexite constante 
0(log(n)) : complexite logarithmique 
O(n) : complexite lineaire 
0(n 2 ) : complexite quadratique 
0(n 3 ) : complexite cubique 
0(2") : complexite exponentielle 
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Complexity 


□ Exemples de temps d’execution en fonction de la 
faille de la donnee et de la complexity de 
I’algorithme. 

On suppose que I’ordinateur utilise peut effectuer 
1 0 6 operations a la seconde (une operation est de 
I’ordre de la |Js) 
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Complexite 






n\T(n) 

log n 

n 

n log n 

n 2 

2 n 

10 

3 \is 

10 |js 

30 |js 

100 \is 

1 000|as 

100 

7 [is 

100 |js 

700 |js 

1/100 s 

10 14 

siecles 

1000 

10 |js 

1 000|is 

1/1 00|is 

1 s 

astrono 

mique 

10000 

13 |js 

1/1 00|is 

1/7 s 

1 ,7 mn 

astrono 

mique 

100000 

1 7 |js 

1/10 s 

2s 
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Complexity 



> Un algorithme est dit polynomial si sa complexity est en 0(n p ). 

^ Un algorithme est dit praticable s’il est polynomial (p < 3). 

Les algorithmes polynomiaux ou p > 3 sont consideres comme tres lents 
(un algorithme polynomial de I’ordre de n 5 prendrait environ 30 ans pour n=1 000) 
^ Un algorithme est dit exponentiel si sa complexity est superieure a tout polynome. 

^ Deux grandes classes de la complexity : 

- IP classe des algorithmes polynomiaux 

- N IP classe des algorithmes « Non determinate polynomiale » 

On a : 

(0(log n) c O(n) e 0(n log n) e 0(n 2 ) c 0(n 3 ) c 0(2") c 0(e n ) e O(nl) 
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Complexity 


Proprietes 

En utilisant la notation de Landau (pour les fonctions de IN dans IR + ), on a : 
f + O(g) = {f + h/he O(g)} 

(h = f + O(g)oh-fe O(g). 
f O(g) = {f h /he O(g)} 
f = Off) 

f = O(g), g = O(h) => f = O(h) 
c Off) = Ofc f) = Off) ( c>0) 

Off) + Ofg) = Off + g) = Ofmaxff, g)) 

Off) + Off) = Off) 

Off) Ofg) = Offg) 
f = Ofg) O g = H (f) 
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Complexity 


□ Calcul de la complexite: regies pratiques 

1. la complexite d’une suite destructions est la 
somme des complexites de chacune d’elles. 

2. Les operations elementaires telle que I’affectation, 
test, acces a un tableau, operations logiques et 
arithmetiques, lecture ou ecriture d’une variable 
simple ... etc, sont en 0(1 ). 

3. T(si C alors A1 sinon A2) = max(T(C),max(T(Al ),T(A2))) 
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Complexity 


4. T(pour i:=e n a e 2 faire Ai fpour) = Y t ( ai ) 

i = e 1 

( si Ai ne contient pas de boucle dependante de i 
et si Ai est de complexity O(m) alors la complexity 
de cette boucle « pour » est 0((e 2 — e n H- 1 )m). ) 

5 . La difficulty, pour la boucle tantque, est de 
determiner le nombre d’iteration Nb_iter (ou 
donner une borne superieure de ce nombre) 

T(tantque C faire A ftantque) = 0(Nb_iter x (T(C) + T(A)) 


SMI Algoll 




Complexite 


□ Exemples 

l. Calcul de la somme 1 +2+...+n 


S: = 0; //0( 1 ) 

Pour i:— 1 an faire 
S:— S + i; // 0(1 ) 


O(n) 


fpour; 


0(1) + O(n) = O(n) 


T(n) = O(n) 
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Complexite 


2. Calcul de : ^w=v / ‘pour i = 2, ... ,n 

j=i 

pour i := 1 an faire 
S := 0; // O(i) 
pour j := 1 a i faire ~ 

S := S + |; // O(i) 

fpour; 

T[i] := s; //O(i) 
fpour; 



T( n ) = 0(n 2 ) 
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3. Analyse de I’algo. Suivant 


Donnee n; (n>0) 
Resultat cpt; 

debut 
cpt := 1 ; 

tantque n > 2 faire 
n := n div 2; 
cpt := cpt +1 ; 
ftantque 

fin 


Que calcule cet algo? 
Quelle est sa complexite? 


Complexite 


□ cpt = 1 + le nombre d’iterations 

□ Le nombre d’iterations = nombre 
de division de n par 2. 

□ Soit p ce nombre. 

- si n est une puissance de 2, i.e. 
n = 2 P alors p=log 2 (n). 

- p verifie: 2 P < n < 2 P+1 

p < log 2 (n)< p+1 =>p=E(log 2 (n)) 

□ cpt = 1 + E(log2(n)), cette 
expression de cpt correspond au 
nombre de bits necessaires pour 
representer I’entier n. 

□ T(n) = 0(log(n)) 
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TRIS ITERATIFS 


SMI ALGOII 1 

V J 



□ Le probleme: 

Etant donnee une suite de n nombres (a lf a 2 ,..., a n ), on 
cherche une permutation (arrangement) des elements 
de cette suite (a’ lf a’ 2 ,..., a’ n ) telle que a’ 1 <a’ 2 ^...^a’ n . 

♦ A partir de la suite (7,1 ,2,6), un algorithme de tri 
donne comme resultat la suite (1, 2,6,7). 

> On se limite aux nombres entiers ranges dans un 
tableau A a n elements. 

> Dans le cas ou les elements sont des collections de 
donnees (enregistrement), on trie le tableau suivant une 
cle (champ de I’enregistrement). 
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□ Tri par selection: 

Alaorithme : 

pour i= 1 a n-1 faire 

-chercher le l eme minimum , A k , de {A ;+7/ ...,A n } 
(K est r indice de m/n{A /+J/ ...,A n } dans le tableau A) 

- echanger A t et A k 

L’algorithme fonctionne selon le schema suivant: 


echange 


A 


trie 





1 i 

SMI ALGOI 

k 


n 




o 
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Analyse du tri par selection 



□ Algorithme: 

Tri_selection(A,n) 

debut 

pour i := 1 a n-1 faire 

/ /Recherche de min{A ; , A n } = A k 

k := i; 

pour /':=/'+ 1 a n faire 

si A[j] < A[k] alors k:=j ; 
fsi; 

fpour 

// ec hange de A k et A t 
temp := A[k ]; 

Mk] := A[i]; 

A[i] := temp ; 

fpour; 

fin 


□ La boucle j determine le i e minimum; 
elle tourne n-i fois (au maximum) pour 
faire (n-i) tests d’elements. 

□ Les echanges de A[k] et A[i] 
demandent 3 operations. 

□ La complexite du corps de la boucle j 
est de la forme 

a(n-i)+b / done en O(n-i) ( i=l , n-1 ). 

□ La complexite de I’algorithme est de 
I’ordre de: 


n-i 

I 

c = l 


(n-i) 


m-l n-l 

i=l i=l 


n(n — 1) 
2 


T(n) = 0(n 2 ) 
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TRI par ISERTION 

□ Algorithme: 

Pour i := 2 a n faire 

- on insere A[i], a sa place, dans le sous-tableau 

A[l..i-1] 

(On cherche le l er element < A[i] parmi {A[i- 1 A[1 ]} en 

decalant d’une position a droite) 

L’algorithme fonctionne selon le schema suivant: 



1 


i-l i 
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Analyse du tri par insertion 



□ Alqorithme : 

Tri_lnsersion( A,n) 
debut 

pour i := 2 a n faire 

cle := A[i]; 
j:= i-1; 

Tantque (j>l ) et (cle<A[j]) faire 
A[j+1]:= A[j]; 

i:= i-i; 

ftantque 

A[j+1 ] = cle; 

fpour 

fin 


□ La boucle j tourne, dans le 
pire des cas, (i-1 ) fois 

(i-1 comparaisons et i-1 
decalages) 

□ La complexite du corps de la 
boucle i est de la forme 

a(i-l)+b, pour i=2,...,n. 

□ La complexite de I’algorithme: 


T g(i -i) + t? 

r=2 


w-l 

^ ak + & 

k=l 


n(n - 1) 

A- T 


+ &(n — 1) 


T(n) = 0(n 2 ) 
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Tri a Bulles 



□ On dit qu’on a une inversion s’il existe (i,j) tels que i<j et 


(ai,...,cii, a i+1 ,...,a n ) 1 =c> 

□ Un tableau est trie s’il n’a aucune inversion. La complexity du 
tri est proportionnelle au nombre d’inversions qui est de 
I’ordre de c n (nombre de couple(i,j) tels que i<j). 


□ L’algorithme consiste a parcourir le tableau a trier en 
examinant si chaque couple d’elements consecutifs (Oj, a i+1 ) 
est dans le bon ordre ou non, si ce couple n’est pas dans le 
bon ordre on echange ses elements et ce processus est 
repete tant qu’il reste des inversions a faire. 
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Tri a Bulles 


□ Algorithmel: 

Bullesl (A,n) 

Debut 

fini := faux; 

Tant que non fini faire 

i: = 1 ; 

tant que (i < n) et (A[i] <A[i+l] ) faire 
i := i + 1 ; 

ftantque ; 
si i < n alors 



echanger (A[i], A[i+ 1 ]) ; 
fini := faux ; } 

s/non 

fini := vrai ; 

fsi; 


ftantque; 



fin 


□ On remarque que les transpositions successives 
font pousser le maximum a la derniere position 
du tableau si on fait un balayage de gauche a 
droite et le minimum a la 1 ere position si on fait 
un balayage de droite a gauche. 

□ Algorithme2: 

Bulles2(A,n) 

Debut i := 1; 

tantque i <n-l faire 

pour j := n a i+1 pas - 1 faire 
si A[j- 1 ] > A[j] alors 
ec hanger(A[j], A[j-1 ]); 
fsi; 

fpour; 
ftantque ; 

Fin 

Le nombre de comparaison d’elements de A < LILli. — 
Le nombre d’echanges < 

T(n) = 0(n 2 ) 
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Exemple 

rw 



Der ech=2 



1 = 2 

Der ech=4 





ALGORITHMIQUE II 


Recurrence et Recursivite 
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Recurrence 


□ Suite recurrente: la definition d’une suite est la donnee 

• d’un terme general defini en fonction du (ou des) terme(s) precedant(s) 

• D’un terme initial qui permet d’initialiser le calcul 

□ Principe de recurrence : 

Soit P un predicat (ou propriete) sur IN qui peut etre soit vrai soit faux (on ecrira 
souvent P(n) a la place de P(n) = vrai). 

On suppose que 

■ P(0) vrai 

■ Vn e IN, P(n) =^> P(n+1) 

Alors , pour tout n e IN, P(n) est vrai. 

Si on considere le predicat suivant 

P(n) : je sais resoudre le probleme pour n 

alors le principe de recurrence nous dit que si je sais resoudre le Pb pour n=0 
et que si je sais exprimer la solution pour n en fonction de la solution pour n+1 
alors je sais resoudre le Pb pour n’importe quel n. 
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□ Examples: 

1 . Puissance 
a 0 = 1 
a n+1 = a a n 

2. Factoriel 

0 ! = 1 

n! = n (n-1 )! , n >1 

3. Suite de Fibonacci 

F 0 = F, = 1 


Ou bien 


n>2 


Recursivite 



i 



a a"' 1 n>0 
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Recursivite 


Un algorithme (ou fonction) est dit recursif s’il est defini en fonction de lui-meme. 
Exemples 

Fonction puissancefx : reel, n : entier) : reel 
debut 

si n = 0 alors retourner 1 
sinon retourner (x * puissancefx , n-1 )) 
fin 

Factoriel (n) 
debut 

si n = 0 alors retournerfl ) 
sinon retourner (n*factoriel(n-l )) 
fin 
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■ fact (n) 
debut 


fin 


si n — 0 alors retourner(l ) 
sinon retourner(n*fact(n- 1 )) 
fsi 


Recursivite 


□ Le deroulement de I’appel de fact(3): 
3*fact( 

3 * 2=6 \ ^ 2*fact(l 


Fact(3) 

6 



*fact(0) 

* 1=1 


La condition n = 0 est appelee test 
d’arret de la recursivite. 


□ II est imperatif de prevoir un test d’arret 
dans une fonction recursive, sinon 
I’execution ne s’arretera jamais. 

□ L’appel recursif est traite comme n’importe 
appel de fonction. 
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Recursivite 


□ L’appel d’une fonction (recursive ou non) se fait dans un 
contexte d’execution propre (pile d’execution), qui 
contient : 

L’adresse memoire de I’instruction qui a appele la fonction 
(adresse de retour) 

Les valeurs des parametres et des variables locales a la 
fonction. 

L’execution d’une fonction recursive se fait par des 
appels successes a la fonction jusqu’a ce que la 
condition d’arret soit verifiee, et a chaque appel, les 
valeurs des parametres et I’adresse de retour sont mis 
(empiles) dans la pile d’execution. 
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Recursivite 


L’ordre des instructions par rapport a un appel recursif est important. 

Exemple: 

afficher(n) 

debut 


si n > 0 alors 

i * 

afficher(n div 10) 


fsi 


r 


<r~ 


ecrirefn mod 10) 


fin 


/.’algorithme recursif afficher(n) permet d’afficher les chiffres d’un entier, strictement positif, selon la 
disposition de I’instruction ecrie(n mod 1 0): 

- Si I’instruction est placee en I , les chiffres sont affiches dans I’ordre inverse 

- Si elle est placee en T, alors les chiffres seront affiches dans le bon ordre 
Pour n = J 23, on a : 

)—> 3 2 1 

r > i23 
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Type de recursivite 


□ Recursivite simple: Une fonction recursive contient un 
seul appel recursif. 

□ Recursivite multiple: une fonction recursive contient 
plus d’un appel recursif (exemple suite de 
Fibonacci). 

□ Recursivite mutuelle( ou croisee): Consiste a ecrire 
des fonctions qui s’appellent I’une I’autre. 

Exemple 
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Recursivite mutuelle 


Pair(n) 

debut 

si n = 0 alors 
retourner vrai 

sinon 

retourner (impair(n-l )) 
fsi 
fin 


Impair(n) 

debut 

si n = 0 alors 

retourner (faux) 

sinon 

retourner (pair(n-l )) 
fsi 
fin 
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Un peu de Structures de Donnees 


□ Notion de pile. 

Une pile est une structure pour representer une suite 
d’elements avec la contrainte qu’on ne peut ajouter, 
ou enlever, un element que d’un meme cote de la 


ajouter 

som 


pile (dit sommet de la pile). 

□ Exemple pile d’assiettes. 

□ Une pile peut etre representee 

par un tableau et un indice de sommet 


C 

nlev 
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Notion de Pile 


□ Operations definies sur les piles: 

initialiser(p : Pile) / /Cree une pile vide. 

sommet(p : Pile) : element// Renvoie I’element au 
sommet de la pile p, sous la condition aue p soit non 
vide. 

empiler(x : element, p : Pile) / / ajoute x au sommet 
de la pile p. 

depiler(p : Pile) / / supprime I’element au sommet de la 
pile p, sous la condition aue p soit non vide . 

pileVide(p : Pile) : booleen / / retourne vrai si p est 
vide. 
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Notion de Pile 


□ Exemple. 

Une expression e est dite bien parenthesee (on se limite au £ (‘ et ‘)’) si : 

1. Le nombre de parentheses ouvrantes (|e| ( ) est egal au nombre de 
parentheses fermantes (|e|j) dans e. 

2. Tout prefixe (partie gauche) u de e verifie: |u| ( - |u| ( > 0. 

Algorithme: on parcourt I’expression e (de gauche a droite). A chaque 

rencontre d’une parenthese ouvrante on I’empile, et a chaque rencontre 
d’une parenthese fermente on depile. 

Si on arrive a la fin de e avec une pile vide, I’expression e est bien 
parenthesee sinon e n’est pas bien parenthesee. 


- I’expression (()())() est bien parenthee. 

- I’expression ())( n’est pas bien parenthesee. 
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Transformation du recursif en iteratif : 
« Derecursivation » 

□ Schema d’un algorithme recursif: 

algoR(X) 
debut 

A 

si C(X) alors 

B; 

algoR(cp(X)) ; 

D; 

sinon 

E; 

fsi; 
fin 


Ou : 

X : liste de parametres 

C : condition d’arret portant sur X 

A, B, D, E : bloc destructions (eventuellement vide) 

<p(X) : transformation des parametres 
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Transformation du recursif en iteratif 
« Derecursivation » 



□ Algorithme iteratif equivalent. 


algoR(X) 


algol(X) 

Debut 


p : Pile 

debut 


A 

initialiser(p); 

si C(X) alors 

A; 


B; 

tantque C(X) faire 
B • 


algoR(<p(X)) ; 

empiler(X, p); 


D; 

X := Cp (X); 

sinon 


A; 

ftantque; 


• ** 

LU 

E; 

fsi; 


tantque (non pileVide(p)) 

fin 


X := sommet(p); 
depiler(p); 


ftantque 

fin 
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Derecusivation 


□ Exemple. 


afficherR(n) 

debut 

si n > 0 olors 
a fficher(n div 1 0) 
ecrire(n mod 10); 
fsi 
fin 


A = B = E = 0 


afficherl(n) 
p : Pile; 

Debut 

initialiser(p); 
tanque n > 0 faire 
empiler(n / p); 
n := n div 1 0; 
ftantque 

tantque (non pileVide(p)) fai 
n := sommet(p); 
depiler(p); 
ecrire(n mod 1 0); 
ftanque 
fin 
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Transformation du recursif en iteratif : 

« Derecursivation » 

□ Recursivite terminale: 

La recursivite est dite terminale si la derniere instruction 
executee est un appel recursive; (Cas ou D = 0). II est 
claire, dans ce cas, d’eliminer la pile dans la version 
iterative. (On depile pour ne rien faire dans la 2 eme 
boucle). 

□ La recursivite d’une fonction F(X) est aussi dite terminale 
lorsqu’elle se termine par I’instruction retourner(F((p(X))). 
On ajoute, dans ce cas, un parametre a la liste X pour 
contenir le resultat de retour d’un appel recursive, 
comme le montre I’exemple suivant: 
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Exemple 



fonction FACR(n) ; 

debut 

si n= 0 alors retourner (1 ) 
sinon retourner (n* FACR(n-l)) ; 
fin. 

RECURSIVE 

fonction FACI(n) ; 
debut 


fonction FACR'(n, r) 

debut 

si n = 0 alors retourner (r) 
sinon retourner (FACR'(n-7, n*r)) 
fin. 

Ny RECURSIVITE TERMINALE 

fonction FACI^n, r) ; 

debut 


r: = 1 ; 

tant que n > 0 faire < 

{ r := n* r ; n := n - 1 ; } 

retourner (r) ; 
fin. 

ITERATIVE 


tant que n > 0 faire 

{ r : = n* r ; n : = n - 1 ;} 

retourner (r) ; 
fin. 
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Complexite des algorithmes recursifs 


complexite des algorithmes recursifs est souvent exprimee 
par une equation de recurrence. 

□ Exemplel. Complexite de I’algorithme recursif pour 
calculer n! (I’operation dominante est la multiplication) 

Soit T(n) le cout de FACR(n). T(n) verifie I’equation: 

T(0) = 0 

T(n) = T(n-l) + 1 

la solution de cette equation est : 

T(n) = n = 1 +1 + ... + 1 (n fois) 
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Complexite des algorithmes recursifs 


□ 

□ 


Exemple2. 

Tri par selection 
sel_rec(T,n) 
debut 

si n > 1 alors 

k<— max {i e{l,2,...,n} / T[i] > T[j] , j=l,2,...,n et j ^ i} 

echanger(T[k],T[n]) 

sel_rec(T,n-l ) 


fsi 


/ 


fin 


□ Complexite : T(n) verifie : 

T(1 ) = 0 

_T(n) = T(n-l) + n 

T(n) = n+T(n-l) = n+(n-l )+T(n-2) =...= n+(n-l )+...+2 +C(1) = 


n(n + 1) 
2 


- 1 


T(n) = 0(n 2 ) 
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On ouvre une parenthese 
Encore un peu de structures de donnees 

□ Notion d’arbre binaire 


On introduit cette notion pour savoir interpreter les arbres d’appels 
dans le cas d’une recursivite double(ou il y a deux appels recursifs). 

- Les arbres sont utilises pour representer une suite d’elements. 

- un arbre est un ensemble de noeuds, chaque noeud represente un 
element de la suite. 


Definition recursive cTun arbre binaire : 

□ un arbre binaire est: 

Soit vide 
Soit forme : 

■ D’un noeud (appele racine) 


arbre de racine r 



■ D’un sous-arbre gauche , note SAG, qui est un arbre binaire) 

■ D’un sous-arbre droit, note SAD, qui est aussi un arbre binaire 
(les deux sou-arbres gauche et droit sont disjoints). 
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Arbre binaire: terminologie 


Soit A un arbre binaire de racine r. 

□ La racine du SAG (resp. SAD) de A 
est appele fils gauche (resp. fils 
droit) de r et r est appele pere. 

Un noeud qui n’ a pas de fils est 
appele feuille. 

un chemin de I’arbre est une suite 
de noeuds n 1/ n 2/ ...n k ou n i+1 est un fils 
(gauche ou droit) de n; , 1 <i<k-l . 

Longueur d’un chemin = nombre de 
noeuds, constituant le chemin, - 1 

Une branche est un chemin de la 
racine a une feuille. 


□ 


□ 


□ 


□ 


□ La hauteur d’un arbre est la longueur 
de la plus longue branche de 
I’arbre. 




-1 si A = 0 


h(A) = 


1 + max(h(SAG(A), h(SAD(A)) 



Branche: 


Branche plus longue 
Hauteur = 4 
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Arbre binaire 


□ Resultat (utile pour la complexity sur les arbres 
binaires de recherche): 


la hauteur h cTun arbre binaire de faille n (n est le nombre de noeuds de I’arbre) 
verifier 1 + flog 2 n] < h < n 

- la hauteur est, en moyenne, un 0(log n) et 
un O(n) dans le pire des cas. 


□ Les algorithmes sur les arbres binaires se ramenent 
souvent aux algorithmes de parcours de ces arbres. 
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Algorithmes de parcours 
puis on ferme la parenthese 


Parcours(a : Arbre) 
debut 

si a ^ 0 alors 
1 


parcours(SAG(a)); 


T 




Traiter(racine(a)) 


parcours(SAD(a)); 

T 


<- 


fsi; 


fin 


Selon la disposition de I’instruction 
traiter(racine(a)), on distingue 3 types de 
parcours: 

) : parcours prefixe. 

T : parcours infixe. 

T : parcours postfixe. 



On passe 3 fois sur un naeud . 

Parcours prefixe: on traite un noeud lorsqu’on le rencontre pour la 1 er fois. 
10-5-9-50-30-20-40-35- 

Parcours infixe: " " " “ la 2e fois. 

5-9-1 0-20-30-35-40-50-60 

Parcours posfixe: " " " le quitte pour la derniere fois. 

9-5-20-35-40-30-60-50- 
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Recursivite double 

& 

Arbre des appels recursifs 


□ Exemplel : calcul de la hauteur 
d’un arbre binaire. 

h(a : arbre) 

debut 

si a = 0 alors retourner -1 
sinon 

hi := h(SAG(a)); 
h2 : h(SAD(a)) ; 
retourner(l +max(h1 ,h2)); 
fsi; 
fin 


□ Un arbre est donne par sa racine 

□ appel de h(l ): 
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Recursivite double 

& 

Arbre des appels recursifs 


Exemple2 : tours de Hanoi (Occupation des moines de HanoY) 

Le jeu consiste a faire passer les 
disques de la tour T1 a la tours 
12 , en ne depla^ant qu’un seul 
disque a la fois, et en utilisant la 
tour intermediate T3 de telle 
sorte qu’a aucun moment un 
disque ne soit empile sur un 
disque de plus petite dimension. 

La solution semble difficile, et pourtant une solution recursive existe. 

Soit n le nombre de disques a deplacer. Si n=l la solution est triviale. 

Si on sait transferer n-1 disques alors on sait en transferer n. 

II suffit de transferer les n-1 disques superieurs de la tours T1 vers la tours T3, 
de deplacer le disque le plus grand de T1 vers T2, puis de transferer les n-1 disques 
de T3 vers T2. Ceci se traduit par I’algorithme recursif suivant: 
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Tours de HanoT 



H(n,n,T2,T3) 

debut 

si n = 1 alors ecrire(Tl, >’,T2) 
sinon 


H(n-1 ,T1 ,T3,T2); 
ecrire(Tl T2) 
H(n-1,T3,T2,T1); 
fsi 
fin 
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H(n,Tl ,T2,T3) 
debut 

si n = 1 alors ecrire(Tl, >’,T2) 
sinon 


H(n-1 ,T1 ,T3,T2); 
ecrire(Tl >’,T2) 


H(n-1,T3,T2,T1); 

Arbre des appels de H(3,a,b,c) 


H (3, a, b, c) 



H (2, a 

, c, b) 

a -> b 

H(2,c 

:,b,a) 







-> b 


H(1 ,a,b,c) 

a — » c 

H(1 ' b ' c ' a ) H(1 ,c,a,b) 

c 

H(1 ,a # b,c) 


a — » b 


b — » c 


c — > a 


a -» b 
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Tours de Hanoi* : Complexity 



Soit T(n) le temps pour deplacer les n disques. T(n) verifie I’equation 


" T(l) = 1 
T(n) = 2 T(n-l) +1 


On a : 

T(2) = 2+1 

T(3) = 2(2 +1) + 1 = 2 2 + 2 + 1 
On montre, par recurrence, que 

T(n) = 1 + 2 + ... + 2 n_1 = 2 n - 1 

Sachant que T(1 0) = 2 10 —1= 1 023 et une annee = 0.3 xl 0 8 secondes, il 
faudrait, pour les moines, 10 10 siecles pour pouvoir deplacer 64 disques! 
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Methode « Diviser pour Resoudre » 


E. CHABBAR 
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DVR 


□ Les algorithmes sont regroupes en families selon 
certains concepts t. q. Division pour Resoudre, 
Glouton, Programmation dynamique. 

□ - L’aspect DVR consiste a diviser le probleme initial 
(de taille n) en sous-problemes similaires de tailles 
plus petites (en general de taille n/2, n/3,...). 

- Ces sous-problemes sont resolus recursivement. 

- On peut combiner ces solutions recursives pour 
avoir la solution du probleme initial. 
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DVR 



Calcul de x n (n > 1 ) 

Le nombre de multiplications par la methode classique ( x n = x x n_1 ) est de I’ordre 
de n. 

La methode DVR exploite la definition suivante: 

1 si n=l 

x n = x p . x p si n= 2p 

x . x p . x p si n = 2p+l 

- pour calculer x n on fait appel (recursif) au calcul x p ou p = E(n/2). (algo. Slide 
suivant) 

- si T(n) est le nombre de multiplications pour calculer x n alors T(n) = T(n/2) + 1 si n 
est pair ou T(n) = T(n/2) + 2 si n est impair. On a, dans tous les cas, T(n) = T(n/2) + 
0 ( 1 ). 

Comme il y a log 2 n divisions euclidiennes successives (avant d’avoir 1 comme 
quotient) , done T(n) = 0(log(n)). (pour s’en convaincre prendre n une puissance de 
2 ) 

(T(n)=0(l )+0(l )+...+0(l ) log 2 n fois) 
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DVR 

Exemples 


i. Calcul de x n (n > 1 ) 

puissBin(x,n) 

Debut 

si n = 1 alors retourner x 
sinon 

si n mod 2 = 0 alors 

y := puissBin(x / n/2) 
retourner (y*y) 

sinon 

y := puissBin(x / n/2) 
retourner (x*y*y) 
fsi; 
fsi; 

Fin. 
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DVR 

Exemples 


2. RECHERCHE DICHOTOMIQUE 

□ La recherche sequentielle d’un element x dans un 
tableau A a n elements est de I’ordre de n. 

□ La recherche dichotomique exiae que le tableau 
soit trie. Elle consiste a: 

- comparer x a I’element du milieu 

- si c’est different, x peut se trouver soit dans la 
moitie gauche soit dans la moitie droite du 
tableau A, selon que x < A[milieu] ou x> A[milieu]. 
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DVR 

Exemples 

2. Recherche dichotomique 

rechDichofAjinfjSupjX) 

debut 

si inf < sup alors 

m := (inf + sup)/2 
si x = A[m] retourner m; 

si x < A[m] alors retourner rechDichofA^n^m-l ,x); 
sinon retourner rechDicho(A,m+l ,sup,x); 
fsi; 
fsi; 

sinon retourner 0; 
fsi; 
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DVR 

Exemples 


2. Recherche dichotomique 

□ Complexite de la recherche dichotomique 

Dans le pire des cas (i.e. x n’est pas dans le tableau) 
le nombre de comparaisons T(n), pour chercher x dans 
un tableau A[1 ..n] a n elements, verifie: 

T(l)= 1 

T(n) = T(n/2) + 1 , n>l 

- la solution de cette equation depend du nombre 
d’iterations (nombre de divisions par 2) 
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DVR 

Exemples 


Iteration 

0 

1 

2 

3 


Nbre d’elts du ss tableau 
n 

n/2 

n/4 

n/8 


p n/2 p 

L’algorithme utilise p iterations(et aussi p comparaisons) et s’arrete 
lorsque n/2 p =1, c.a d. p = log 2 n. 

par consequent, la recherche dichotomique est en 0(log 2 n), i.e. 

T(n) = 0(log 2 n) 
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3. Tri par fusion 

■ Fusion de deux tableaux tries. 

Etant donnes deux tableaux tries T1 [1 ..nl ] et 
T2[l ..n2]. La fusion consiste a construire un tableau 
T[1 ..nl +n2] contenant tous les elements de T1 et 
T2 dans I’ordre croissant. 


T1 T2 




Fusion 



Fusion(Tl ,nl ,T2,n2) 

// T est un tableau qui contient le resultat de la fusion 

i 1 1 ; i 2 : — 1; k:=l; 

tantque (il <nl ) et (i2<n2) faire 
si T1 [il ] < T2[i2] alors 
T[k] T1 [il ]; 
k: = k+l ; il il +1 ; 

sinon 

T[k] := T2[i2]; 
k: — k+ 1 ; i2 := i2 + 1; 
fsi; 

ftantque; 

// on recopie les elts restants dans I’un des //tableaux Ti dans le tableau T 

si il < nl alors copier(Tl ,il ,nl ,T,k) 
sinon copier(T2,i2 / n2,T / k) 
fsi; 

retourner (T); 


copierfAjd^Bjj) 

debut 

pour i:— d a f faire 

m S= A[i] 

i*= i+1; 

fpour; 

retourner(B); 

fin 
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la complexite de la fusion est en Ofnl +n2) 




TriFusion 


□ Le tri par fusion est un exemple typique de la 
strategie diviser pour resoudre, a savoir: 

1. Diviser le tableau T[1 ..n] en deux sous-tableaux 
T[l..E(n/2)] et T[E(n/2)+l ..n] 

2. Trier (recursivement) les deux sous-tableaux (2 
appels recursifs a la meme fonction TriFusion). 

3. Fusionner les deux sous-tableaux; 

Ceci est schematise par I’exemple suivant: 
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TriFusion 


Diviser 


Resoudre recursivement 


Combiner 
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TriFusion 



TriFusion(T, n) 

// T1 , T2 : des tableaux (locaux) de longueurs variables a chaque appel 

debut 


si n > 1 alors 

copier(T, 1 , n/2, T1 ,1 ); 
copier(T, n/2 + 1 , n, T2 , 1 ); 

T1 := TriFusion(Tl , n/2); 

T2 := TriFusion(T2, n — n/2); 

T := Fusion(Tl, n/2, T2 , n - n/2); 


fsi; 

retourner(T); 



fin 


Retourner(Fusion(TriFusion(T 1 ,n/2),n/ 2,TriFusion(T2 / n-n /2) / n-n / 2) 
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TriFusion: exemple dexecution 

Tableau a trier: 


3 

9 

1 

6 

5 

2 

4 

8 



TF(1 ..4) 


3 

9 

1 

6 


TF(1 ..8) 



1 ..4,5. .8) 


8 



TF(l..l) TF(2..2) F( 1 -.2) 


TF(1 

~2) 


TF(3 


FM 

..2,3 

-.4) 

3 

9 


1 

6 

1 

3 

6 

9 


TF(3..3) TF(4..4) F(3..4) 


1 


TF(5..6) 


TF(7..8) 


8 


TF(5..5) TF(6..6) F(5..6) TF(7..7) TF(8..8) 



i \ . 

2 

4 

5 

8 

) F(7..8) 


4 

8 



Recopie de ss-tableau 



Tableau (vert)retourne 
a I’appel recursif 


Fusion des 2 ss-tableaux(rouges) 
resultats des appels recursifs 
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DVR: Complexite 

□ La complexite T(n) pour trier un tableau de n 
elements par I’algorithme TriFusion verifie: 

; T(i ) = o 

T(n) = 2 T(n/2) + O(n) , n > 1 


(II y a 2 appels recursifs, chacun porte sur la moitie du 
tableau. O(n) pour recopier les 2 ss-tableaus en 
2xO(n/2) + leur fusion en O(n))). 
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Equation de recurrence des alg.DVR 


□ La recurrence, utilisee par les algorithmes type DVR, est 
souvent de la forme: 

0(1 ) pour n=l 



a T(n/b) + 0(n d ) , n > 1 


a : est le nombre de divisions du probleme initial en sous-probleme 
(nombre d’appels recursifs) 


n/b : la taille de chaque sous-probleme (b>2) 


0(n d ) : le temps necessaire pour decomposer le probleme en sous- 
probleme + le temps pour combiner les solutions des ss-problemes 
pour avoir la solution du probleme initial. 
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DVR: Resultat de Complexity 


□ Theoreme: 

Soit T : IN IR+ une fonction croissante a partir d’un certain rang, telle 
existe des entiers n 0 >l, b>2 et des reels d>0, a>0 pour lesquels 


T(n 0 } - fe 

T(n) = aT + cn* 


Alors on a : 


T(n) = 


O (n d ) 


0(71 d l0g t 7l ) 

O(n logba ) 


si a 


71 

n > — puissance de b 

7Lq 


< b d 


SI 


i a = b d 


si a > b d 


SMI-ALGOII 




DVR: Resultat de Complexity 


sort T{n\ = ar(^) + . 

On montre par recurrence sur p (p>l ), que : r OO = + cnd ^ (^V 


y-i 


i=0 


prenons — = b p f et done p = 0(log b n ).Laretationpredente devtent : 
n Q 

L p - i 

T( it) = k& p -+■ cn d ^ (^) J ' = f(ri) + cn d ^(n) , oi : 

i-o 

f (n) = *aP = tf(n lo “ fl ) et g{n) = 


dlli/altfruffi a P = e ?lna = e log& Hc lnC; = ff lx ^In b =(— )l°E6 a = 0(n log:£!a ) 

a f"P— i 

g{n)est une suite geometrique de raison r = -T J par suite g(n) = (r 1) 


Irr 


In t 


n 


Z 7 d 


r-1 


premier cas: sir= = 1, alors g{n} = p= O(log & n) et par consequent = 0[n d log£?l} si 6 = b d 

dmixismecas :g{n) = — < sir = 7-7 < 1 etg^n) = 0(1) dans ce cas , p&r consequent rOO = si a < & £ 

1 — r 1 — r b a 


^ 

traisihne cas-.g(n) = ■< r p sir = . d >1 etg(y%) = o(n~ d r^ OSba ^du faite que: 


a 


b d 


fTou T[n) = o( n l ° Ej a } si a > b' 


a? 

r'P — 

b* d 
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"Jtl 

= 0(y°-s& E )(— )~ d = 0(n lo ^ a )0(n- d ) 

n a 




DVR: Resultat de Complexity 


D’une maniere generale. Si 

T(1 ) = 0(1 ) 

T(n) = a T(^ ) + 0( n rf ) ( n > 1 ,a> 0 , b> 1 , d > 0) 

Alors 


0(n d ) 


si a < b d 



0(n d log b rc ) 




si a = b d 
si a > b d 
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DVR 


□ Applications 

1. Pour la recherche dichotomique et le calcul de la 
puissance on a : T(n) = T(n/2) + 0(1 ) 

done T(n) — 0(log n) (a= 1 , b=2, d=0) 

2. Pour le TriFusion on a : T(n) = 2 T(n/2) + O(n) 

Alors T(n) = 0(n log 2 n) (a=2, b=2, d= 1 ) 

3 . T(l) = 1 

T(n) = 2t(n/2) + 0(n 2 ) 
a pour solution T(n) — 0(n 2 ) (a=b=d=2) 
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DVR : tri rapide (quicksort) 


□ Soit a trier le tableau T[g..d] (au depart g— 1 et d— n) 

□ Le principe de I’algorithme reside dans une procedure, 
appelee partition, qui reorganise les elements de T 
autour d’un pivot (element du tableau choisi au hasard) 
de sorte que : 

1 ) II existe un indice p (g< p <d) tel que p est la 
position definitive du pivot (T[p] — pivot). 

2) tous les elements T[g], ..., T[p-1] sont inferieurs ou 
egaux a T[p]. 

3) tous les elements T[p+1], ...T[d] sont superieurs ou 
egaux a T[p]. 
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DVR : tri rapide 



□ Le travail de la fonction partition consiste a: 

o Choisir un element du tableau comme pivot(par exemple le premier 

T[g]) 

Parcourir le tableau depuis la gauche(de gauche a droite) jusqu’a 
rencontrer un element > T[g] 

Parcourir le tableau depuis la droite (de droite a gauche) jusqu’a 
rencontrer un element < T[g] 

o Echanger ces deux elements dans le tableau 

o Continuer ce processus jusqu’a ce que les deux indices (de gauche et 
de droite) se croisent. 


o 


Finalement, e changer le piv ot T[g] et I’element indiqu e par I’in dice 
de droite. 


< pivot 


Non encore analyse 


> pivot 
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► i 


I 


g 


d 




DVR : Tri Rapide 


/ /T[n+1 ] = + oc 
Partition(T, g, d) 
debut 

pivot := T[g]; 

i g; i := d+1; 

tantque i < j faire 

i := i +1; 

tantque T[i] < pivot faire i:=i+l; 
ftantque 

i := i -1; 

tantque T[j] > pivot faire j:=j-l; 
ftantque 

si i < j alors echanger(T, i, j); fsi 
ftantque 
echanger(T, g, j); 
retourner(j) 
fin 


□ Exemple 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (1 0) (1 1 ) (1 2) (1 3) 
3141592653 5 8 9 ij 

3141592653 5 8 9 3 10 

3131592654 5 8 9 3 10 

3131592654 5 8 9 57 

3131295654 5 8 9 57 

3131295654 5 8 9 65 

2131395654 5 8 9 5 


1 3 1 


9 5 6 5 4 


8 
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DVR : Tri Rapide 

□ La fonction partition, appliquee a un tableau T, 
produit trois sous-tableaux: 

o Un sous-tableau reduit a un seul element T[p] qui 
garde sa place definitive dans le tri de T, et 

o Deux sous-tableaux T[g .. p-1], T[p+1 .. d]. 

□ Pour trier T, il suffit d’appliquer recursivement le 
meme algorithme sur les deux sous-tableaux. 
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DVR : Tri Rapide 

quickSort(T, inf, sup) 
debut 

si inf < sup alors 

p := partition(T, inf, sup); 
quickSort(T, inf, p-1 ); 
quickSort(T, p+1, sup); 
fsi 
fin 
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Complexite du Tri rapide 


La complexite de la fonction partition, appliquee a 
T[1 ..n], est en O(n). 

- Cas le plus defavorable : 

Cas ou le pivot sort, a chaque fois, en premier 
(ou en dernier) element (T: tableau trie). 

La partition coupe le tableau en un morceau de un 
element et un morceau de n-1 elements, dans ce 
cas on a : 

C(n) = C(n-l) + O(n) 

(O(n) est le coOt de la partition) 

On en deduit que C(n) est en 0(n 2 ) 
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Complexite du Tri rapide 

- Cas le plus favorable : 

cas ou le pivot est I’element mediane de T. 

La partition coupe T en deux morceaux de 
taille n/2 

C(n) = 2 C(n/2) + O(n) 

ce qui donne: C(n) = 0(n log n) 

La complexite moyenne est aussi de I’ordre de n log n 
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Complexite du Tri rapide 



- Complexite moyenne du tri rapide 

La formule de recurrence donnant le nombre de 
comparaisons effectuees par le tri rapide pour une 
permutation aleatoire de n elements verifie 

C 0 = C, = 0 et 


i 


n-l 


C n = "-1 + - ZlCi + C.,,) 

i=0 


pour n >2 

Le terme generique C n peut s’ecrire : 


o n-l 

C n = " “ 1 + TP I Ci 

i=0 
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•Complexite moyenne du tri rapide 


n-1 


• C n n 1 + n iio C ' 


n-2 


C n = n — 1 + — C , + ^ 

n n i=0 


I c, 


C , + 

I 


n-2 


C — n — 1 + v “'n-l 

n n 


n-1 o 

(n-2+ 2 IC,) 

n-1 i=0 


n 


( n-1) (n-2) 
n 


C n = 


c = 


n 


C n + 

n-1 


n-1 


n 


C i + 

n-1 


2n-2 


n 


n+1 


n 


C i + 

n-1 


2(n-l ) 


n 


En posant D n = 


n+1 


On aura la recurrence : 


D = D , + 

n n-1 


n+1 n(n+1) 


En negligeant le dernier terme de D n on a: D n — log n 

Du fat que : 1 + ~7~ + ••• + 1 

n 


— - L Dll 1 2L& < rl Cl C " est en 0(n log n) 




Preuve de programmes 


Notions de Logique 


E. CHABBAR 
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Logique 


La logique a joue un role decisif dans le 
developpement de I’informatique, notamment en 
informatique theorique: 

- Definition d’un modele theorique des premiers 
ordinateurs (Machine de Turing) 

- Calcul booleen pour la conception et I’etude des 
circuits. 

- La recursivite pour definir la calculabilite. 

- la decidability et la complexity pour etudier la limite 
de la machine. 

- La programmation fonctionnelle. 
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Logique 


□ Actuellement I’informatique a envahi tous les 
domaines de la vie. Le probleme le plus important 
qui se pose dans la conception d’une application 
informatique est de prouver qu’elle est exempte 
d’erreurs et qu’elle resout le probleme pour 
laquelle elle a ete congue. Pour cela, on definit une 
tache par une formule de logique et on montre 
qu’elle est satisfaite par un modele de cette 
application. Ce type de preuve est appele « 
Verification Formelle ». 
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Logique 


□ Une logique, par definition, est un ensemble de formules. 

□ Une formule est construite, sur un alphabet, suivant certaines 
regies (Syntaxe). 

□ La semantique d’une formule est une valuation (ou 
interpretation dans un modele) qui determine la valeur de 
verite de la formule. 

□ II y a plusieurs types de logiques: 

- Logique des propositions (d’ordre 0) 

- Logique du premier ordre (les predicats en font partie) 

- logique du second ordre et logique d’ordre superieur. 
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Logique des propositions 


□ - On note P = {p, q, ...} I’ensemble des propositions 
atomiques. Chaque proposition atomique est une 
variable qui ne peut prendre que « vrai » ou « faux ». 

{ 1 / a, •£=>} I’ensemble des connecteurs (ou 

operateurs) logique plus les parentheses. 

□ - On note F(P) I’ensemble des formules deduit de P. 

□ F(P) est le plus petit ensemble qui contient P et qui est 
stable par les connecteurs. 

En d’autres termes: 

Une formule est une suite de symboles de P U {-i , a, 
v, =>, <=>, (, )} construite selon les regies suivantes: 
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Logique des propositions 


- Toute proposition atomique est une formule. 

- si f est une formule alors —if est aussi une formul 

- si f 1 et f2 sont des formules alors : 

- fl a f2 est une formule 

- fl v f2 « « « 

- fl => f2 

- fl <=> f2 « (( (( 

Exemple: (— ip => q) v (p a q) est une formule. 
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Logique des propositions 


□ La semantique des operateurs logiques est donnee 
par des tables de verite. 
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P red i cats 


□ Les predicats sont construits avec : 

- les constantes (0,1,2,....) 

- les variables (x, y, ....) 

- les fonctions (f, g, +, *, ...) 

- des connecteurs logiques (— i , A, v, =>, 

- des parentheses 

- des quantificateurs (V, 3) 

□ Un predicat atomique est un predicat qui ne contient ni 
connecteur ni quantificateur. 

- Exemples: x < y ; pair(2 x) ; 
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Predicats 


Regies de formations des formules de predicats 

- Tout predicat atomique est une formule 

- si f 1 et f2 sont des formules alors 

— i f 1 , f 1 v f2, f 1 Af2,fl =>f2,fl<=>f2 
sont des formules. 

- si f est une formule alors 

V x f est une formule 
3 x f est une formule 
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Predicats 


□ Exemples: 

1 ) Vx pair(x + x) 

2) 3x premier(x) a x < succ(succ(0)) 

3) Vx (oiseau(x) => vole(x)) (tous les oiseaux volent) 

4) 3x (oiseau(x) a — i vole(x)) (4 = — i 3) 

5) Vi Vj (1 <i <j<n — ^ T[i] ^ T [ j] ) (specification d’un tableau 

□ Remarques 

-i (Vx P) = 3x — iP 
-i (3x P) = Vx — iP 
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Deduction logique 


□ Demonstrations logiques 

□ Un sequent est un couple de la forme ($ ,f), ou f 
est une formule et $ un ensemble fini de 
formules . L’ensemble $ est I’ensemble 

des premisses du sequent, la formule f sa 

conclusion. 


SMI-ALGOII 




• Sequents prouvables 

• Un sequent ,f) est prouvable , ce que Ton notera 
$ \- f, s’il peut etre construit en utilisant un nombre fini 


de fois les 6 regies suivantes 


1 . 

2 . 

3 . 

4 . 

5 . 

6 . 


si fe^, alors $ |-f 

sig |-f, alors g|-f 

si S' b ( f => f ’) et S' f , alors $ \- f ’ 

si ^ , f |-f’, alors 

# b f ssi ^ b— f 

si $ ,f b f’ b alors $ h f 


(Hypothese) 


(Modus Ponens) 
(Synthese) 


(Raisonnement par I’absurc 
(Contradiction) 
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Demonstrations logiques 

Une demonstration d’un sequent prouvable $ |- f est 
une suite finie de sequents prouvables f\ |- fi, 

i = 1 ,...,n, telle que : 
jn = j et f n = f 

chaque sequent de la suite est obtenu a partir des 
sequents qui le precede en appliquant I’une des 6 regie 

Remarque : le premier sequent de la suite est 
necessairement obtenu par utilisation d’une hypothese 
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• Demonstrations logiques - Exemple 

• Tous les hommes sont mortels, et 

• les Grecs sont des hommes, done les 
Grecs sont mortels 

h = « etre un homme » 

m = « etre mortel » 

g = « etre Grec >> 

(h m), (g => h) |- (g => m) 
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• Demonstrations Logiques - Exemple 

• Tous les hommes sont mortels, et 

• les Grecs sont des hommes, done les Grecs 
sont mortels 


1 . 

(h => m), (g => h), g 

-g 

hypothes 

2. 

(h => m), (g => h), g 

- (g => h) 

hypothes 

3. 

(h => m), (g => h), g 

-h 

MP (1&2) 

4. 

(h => m), (g => h), g 

- (h => m) 

hypothes 

5. 

(h => m), (g => h), g 

- m 

MP (3&4) 

6. 

(h => m), (g => h) |-(g => m) 

synthese 
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Preuve de programmes 




Logique de Hoare 
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Correction de programmes / 

Specification 


□ un programme est correct s’il effectue la 
tache qui lui est confiee dans tous les cas 
permis possibles 


□ necessity de disposer d’un langage de 
specification permettant de decrire 
formellement la tache confiee a un 
programme. 
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Correction de programmes / 

Specification 

□ description des proprietes que doit 
satisfaire un programme pour repondre au 
probleme pose 

- relation entre les entrees et les sorties 
du programme 


Entree E — 
Precondition : P(E) 


Programme 


— * Sortie S 

Postcondition : Q(E, S) 
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Correction de programmes / 

Specification 


□ Exemple 

une specification pour le probleme du 
calcul de la racine carree entiere par 
defaut : 

Donnees n : entier ; 

Resultat r : entier ; 


Pre condition : n > 0 

Algorithme: r := 0; tantque (n >(r +1) 2 ) faire r := r + 1 ftantque 
Post condition : (r 2 < n) a (n < (r+1) 2 ) 
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Test vs Verification 



• Les couples (n, r) de I’ensemble suivant {(0, 0), (1, 1), (2, ), (3, 1), (4, 2), ...} 
sont des tests qui reussissent 

• Hoare propose une execution de I’algorithme sur une valeur 
symbolique qui est un ensemble de valeurs defini par une expression 
logique (ou predicat). 

• L’ensemble {(0, 0), (1, 1), (2, 1), (3, 1), (4, 2), ...} est 

defini par {(n,r) / n>0 a r >0 a r 2 < n a (n < (r+1) 2 } 

• On note {n>0 a r >0 a r 2 < n a (n < (r+1) 2 } la valeur symbolique des 
variables n et r. 

• L’execution d’un algorithme A sur une valeur symbolique de donnees 
definies par I’expression { p} qui donne une valeur symbolique resultat 
definie par I’expression { q} est notee {p} A {q} . {p} A {q} est appele triplet 
de Hoare (p : precondition, q : postcondition). 
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Systeme formel 


• Interpretation d’un triplet de Hoare 

{p} A {q} signifie : 

Si la propriety p est vraie avant I’execution de A ET si I’execution de A 
se termine, ALORS la propriety q est vraie apres I’execution de A. 
(Correction partielle) 

• Un systeme formel est un triplet <L, Ax, R> ou : 

- L est un langage definissant un ensemble de formules, 

- Ax est un sous-ensemble de L ; chaque formule de Ax est appelee axiome, 

R est un ensemble de regies de deduction de formules a partir d’autres formules 
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Langage algorithmique 

Instructions : 

- Affectation (:= symbole d’affectation et = pour la comparaison) 

- Controle : 

Si <cond> alors instruction fsi 
ou Si <cond> alors instruction 
sinon instruction fsi 

- Boucle : 

Tantque <cond> faire instruction ftantque 
Restriction : 

- Pas de fonction (ou procedure) 

- ni de variable pointeur 

- pas de designation de la forme t[i] ou i fait reference a un autre 
tableau. 
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Logique de Hoare 

• Definition de la logique de Hoare : La logique de 
HOARE est un triplet <L, Ax, R> avec : 

- L est I’ensemble des formules { p} A { q} ou p 
et q sont des predicats et A est un fragment 
d’algorithme(ou de programme) 

- Ax est I’ensemble des axiomes de la forme 
suivante : { p[ x / e] } x := e { p(x)} 

(p(x) est obtenu de p[x/e] en substituant toute 
occurrence de e par x) 

- R est I’ensemble de regies de deduction 
suivant : 
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Logique de hoare 


(Seq) si {p} A1 {r} ; {r} A2 {q} alors 

{p} A1 ; A2 {q} 

(cons g ) si p p’, {p’} A {q} alors 

{P} A {q} 

(cons d ) si {p} A {p’}, p’ q alors 

{p} A {q} 

(condj si {p A c} a {q} , (p A — ic) =^> q alors 

{p} si c alors A fsi {q} 

(cond 2 ) si {p A c } A1 {q},{p A — ic} A2 {q} alors 

{p} si c alors A1 sinon A2 fsi {q} 
(tantque) si {I A C} A {1} alors 

{ 1} tantque C faire A ftantque {I A — iC} 


(I doit etre un invariant de la boucle) 
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Exemple : regies d’affectation 

et consequence 

• { p[ x / e] } x := e { p(x)} 

1. {n+1 > 0} n:= n+1 {n>0} (x=n, e = n+1) 

2. {k > 0} n :=0 {k > n} (x=n, e =0) 

3. {x = 4} x := x+1 {x=5} en ecrivant: 

{(x+1 ) -1 =4} x := x+1 { x - 1 =4} => {x = 5} 

4. {x > 0} x := x + 1 {x > 1 } 

{x > 0} => {x + 1 > 1 } 

X := X + 1 
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Exemple : regies de condition 

5. [vrai} 

si x > 0 alors 

{vrai a x > 0} => {x > 0} =>{x+ 1 > 1 } 

X := X + 1 ; 

{x>l} 

sinon 

{vrai a x < 0} => {x < 0} => {-x > 0} 
x := - x; 

{x > 0} {x > 1 } 
fsi 

{X > 1} 
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Schema de preuve d’un algorithme iteratif 

Tant que 

□ Soit I’algorithme suivant: 

debut 

{pre} 

init; 

{ 1 } 

tantque C faire 
{IaC} 

A 

{ 1 } 

ftantque 

{I A — i C} => 

{post} 
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Exemple: calcul de la racine carree 

par defaut 

■ Algorithme A: 
donnee : n : entier; 
resultat : r : entier; 
debut 
r := 0; 

tantque n>(r+ l) 2 faire 
r := r + 1 ; 

ftantque 

fin 


- Specification: 
precondition : {n > 0} 
postcondition : { (r 2 < n a n < (r + 1 ) 2 ) } 

■ invariant de la boucle: I = {r 2 < n} (evident, sinon on le deduit de la postcondition) 
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Algorithme A annote et prouve 


Algorithme A: 
donnee : n : entier; 
resulted : r : entier; 

Debut 

{n > 0} => {n > 0*0} 
r := 0; 

{n > r*r} 

tantque n > (r + 1 ) 2 faire 

{n > r 2 A n > (r + 1 ) 2 } 
r := r + 1 ; 

{n > r 2 } 

ftantque 

{r 2 < n A n < (r + 1 ) 2 } 
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Exemple: 



■ Algorithme B: 
donnee : n : entier 
resultat : y : entier 
debut 

X := n; 

Y := 1; 

tantque x > 1 faire 

y := y * x; 

X := X — 1 ; 

ftantque 

fin 


■ Specification: 
precondition : {n > 0} 
postcondition : {y = n!} 

■ Invariant: I = { n! = y x! a x > 0} 


calcul de n! 
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Algorithme B annote et prouve 



Algorithme B: 
donnee: n : entier; 
resulted: y : entier; 

debut 

{ n - 0 } 

X := n; 

{x > 0 A x = n} 

y := 1; 

{n! = y x! A x > 0} 

tantque x > 1 faire 

{n! = y x! A x > 0 A x > 1 } => {n! = (y x)(x-l )! A x > 1 } 

y := y * x; 

{n! = y (x-1 )! A x-1 > 0} 
x := x — 1 ; 

{n! = y x! A x > 0} => {n! = y x! A x > 0} 

ftantque 

{n! - y x! A x > 0 A x < 1 } => { y = n!} 

Exercice: Faites la preuve de I’algo. de n! en faisant une boucle qui va de 1 an 
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